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C^ ; ABSTRACT 

Context. Improving our understanding of the complex star formation process in clusters requires studies of star-forming clouds to 
yf^ ' search for dependencies of the physical properties on environmental variables, such as overall density, stellar crowding and feedback 

from massive (proto-)stars. 

Aims. We aim to map at small spatial scales the temperature and the velocity field in the protocluster associated with IRAS 
05345+3157, which contains both intermediate-/high-mass protostellar candidates and starless condensations, and is thus an ex- 
cellent location to investigate the role of massive protostars on protocluster evolution. 

Methods. We observed the ammonia (1,1) and (2,2) inversion transitions with the VLA. Ammonia is the best thermometer for dense 
and cold gas, and the observed transitions have critical densities able to trace the kinematics of the intracluster gaseous medium. 

(-H . Results. The ammonia emission is extended and distributed in two filamentary structures. The starless condensations are colder than 

Oj' the star-forming cores, but the gas temperature across the whole protocluster is higher (by a factor of 1.3-1.5) than that measured 

I I typically in both infrared dark clouds and low-mass protoclusters. The non-thermal contribution to the observed line broadening is 

O ■ at least a factor of 2 larger than the expected thermal broadening even in starless condensations, contrary to the close-to-thermal 

line widths measured in low-mass quiescent dense cores. The NH3/N2H^ abundance ratio is greatly enhanced (a factor of 10) in the 

C/3 ' pre-stellar core candidates, probably due to freeze-out of most molecular species heavier than He. 

Cu ' Conclusions. The more massive and evolved objects likely play a dominant role in the physical properties and kinematics of the pro- 

tocluster The high level of turbulence and the fact that the measured core masses are larger than the expected thermal Jeans masses 
indicate that turbulence likely was an important factor in the initial fragmentation of the parental clump and can provide support 

K^ against further fragmentation of the cores. 

CN ' Key words. Stars: formation - Radio lines: ISM - ISM: individual (IRAS 05345-1-3157) - ISM: kinematics and dynamics - ISM: 

molecules 
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Cf^ 1. Introduction al. l2007h . Such interaction is expected to be very important in 
^^ clusters containing intermediate- and high-mass protostars or 
(N| Most stars form in clusters in the densest portions of giant newly formed massive stars, given the typical energetic feedback 
^H , molecular clouds, but the initial conditions of the star forma- provided by these objects (powerful outflows, strong winds, UV 
L; . tion process in clusters are still poorly understood. Studies have radiation, expanding Hii regions). There is also vigorous theoret- 
unveiled the physical and chemical properties of isolated low- j^al debate on how star formation proceeds in clustered regions, 
mass pre-steUar cores, i.e. stariess condensations on the verge Specifically: is the local kinematics of the gas dominated by 
H of gravitational coUapse, in nearby low-mass star-forming re- feedback from protostellar outflows of already-forming, gener- 
C? . gions. These works show that pre-stellar cores have dense (n ~ ^Uy low-mass stars (Nakamura & Li 5oo7) or by feedback from 
10^ - 10 cm- ) and cold (T ~ 10 K) nuclei, and that their high-mass stars? To what extent does turbulence regulate the star 
internal motions are thermally dominated, as demonstrated by formation rate? Can a model of star formation starting from qui- 
their close-to-thermal line widths, even when observed at low escent starless cores be applied to clustered regions (McKee & 
angular resolution (see Bergin & Tafalla l2007| for a review), jan 2003)? To put constraints on theoretical models it is crucial 
One as yet pooriy investigated major issue is if, and how, the ^^ determine the kinematics and the physical properties of star- 
environment influences the physical and chemical properties of jegg cores in clustered environments, especially those containing 
these pre-stellar cores. In clusters containing several forming sources of very energetic feedback 
stars in a small area (< 0.1 pc), turbulence, relative motions, 

and interactions with nearby forming (proto-)stars can in prin- An observational eff^ort in fliis direction has started, but it is 

ciple affect the less evolved condensations (Ward-Thompson et mosdy concentrated on nearby low-mass star-forming regions 

like Ophiuchus (see e.g. Andre et al. 2007; Friesen et al. 2003), 
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al. (l2009l) show that cores within low-mass star-forming clusters 
have typically higher kinetic temperatures (~ 15 K) than low- 
mass isolated ones (10 K). On the other hand, the kinematics 
seems to be dominated by thermal motions like in more isolated 
cores , even though the exter nal env ironment is turbulent (Andre 
et al. l2007l Rathborne et al l2008l) . Few works have been per- 
formed towards protoclusters forming intermediate- and high- 
mass stars to dat e (Palau et al. 2007a, 2007b, Leurini et al. 2007, 
Pillai et al. l2011h . due to the fact that the typically large distances 
(> 0.5-1 kpc) of high-mass star-forming regions make the study 
of clustered environments challenging. 

The target of the present paper is a well-studied protoclus- 
ter embedded in a massive (180 Mq) pc-scale clump (Fontani 
et al. 120061 Klein et al. |2005|) . The protocluster harbours 
intermediate- to high-mass protostellar objects and starless con- 
densations (Fontani et al. 2008, Fontani et al. 2009, hereafter 
paper I and paper II, respectively), and is located at a distance 
of 1.8 kpc, ~ 1' north-east from the IRAS source 05345+3157 
associated with a rich cluster of infrared sources (Varricatt et 
al. l2010h . Hereafter, the protocluster will be called 105345. One 
can summarise the main findings of the previous observations as 
follows (as general reference, see Fig. [T}: 

- at a resolution of ~ 3", the millimetre continuum reveals 
three dusty cores, CI, C2 and C3. CI is then resolved into 
two peaks, CI -a and Cl-b, down to 1'.'5 resolution (paper I, 
paper II); 

- both CI -a and Cl-b are associated with infrared emission 
and are believed to harbour intermediate- to high-mass pro- 
tostars, with Cl-b being the most massive and associated 
with a hot core, unlike CI -a; in paper II we wrongly claimed 
that Cl-b is likely associated with 3.6 cm-continuum emis- 
sion (Molinari et al. 2002) due to an error in the coordinates 
of the peak. The radio continuum emission is instead more 
likely associated with C2, and it can be due either to a ra- 
dio jet or to a deeply embedded Hii region. In either case, C2 
harbours a very young intermediate- or high-mass object; C3 
is also associated with both dense gas and an infra red sou rce, 
but it s natu re is less clear (paper I, Varricatt et al. 1201 0[ Lee 
etal. l20TTl) : 

- an extended gaseous emission detected in N2H^ (1-0) with 
the PdBI encompasses all continuum cores and reveal addi- 
tional condensations (paper I, paper II); 

- two of these condensations, called N and S, were also de- 
tected in emission lines of N2D^ but not in the mm- and 
IR-continuum. Based on the chemistry and on the lack of 
infrared emission, N and S are probably starless and in the 
pre-stellar phase (paper I); 

- a widespread CO outflow hkely driven by Cl-b and/or CI -a 
interacts with the starless condensations (paper II). 

These features make this protocluster an excellent location to in- 
vestigate if and how starless cores are influenced by the feedback 
from forming intermediate- to high-mass stars. In this paper we 
present new VLA high-angular resolution observations of the 
(1,1) and (2,2) inversion transitions of NH3 towards 105345, fo- 
cusing mainly on the kinetic temperature and the kinematics of 
the dense gas. NH3 is the ideal thermometer because it is known 
to be a tracer of relatively dense gas and to be present in the gas 
phase even in regions where CO and other C-bearing species are 
frozen onto dust grains, and the NH3 (2,2) to (1,1) line ratio is 
sensitive to temperature. Also, the critical density of both transi- 
tions (about lO'* cm"-') allows us to investigate the kinematics of 
the gaseous envelope in which the dense protocluster cores are 
embedded. 



In Sect. |2] we outline the details of the VLA observations 
and data reduction. Sect. |3]presents the immediate observational 
results and the method used to derive the physical parameters 
from NH3, while in Sect. |4] we show the maps of these parame- 
ters. In Sect.|5]we discuss some aspects of the chemistry and the 
implications of this work in the context of the protocluster for- 
mation problem. A summary and the main conclusions are given 
in Sect. |6] 



2. Observations and data reduction 

The NRAO Very Large Array (VLA)0] ammonia observations to- 
wards 105345 were performed on October 31 and November 2, 
2009. As phase centre we used the nominal position of the sub- 
millimetre peak detected with the JCMT (Fontani et al. 2006), 
namely RA(J2000) = 05*'37™52^.4 and Dec(J2000) = 32°00'06", 
which roughly corresponds to the position of CI . The local stan- 
dard of rest velocity Vlsr is -18.4 km s"'. The NH3 {J,K) - 
(1,1) and (2,2) inversion transitions at 23.694496 and 23.722633 
GHz, respectively, were observed simultaneously, using the 2- 
IF spectral line mode of the correlator, with 3.125 MHz band- 
widths: the first IF covered the main line and the two innermost 
satellites of the (1,1) transition, while the other IF was used to 
observe the (2,2) line. Channel spacing was 24.414 kHz, corre- 
spoding to a spectral resolution of ~ 0.3 km s"' at the frequency 
of the lines. The array was used in the most compact configu- 
ration (D), offering baselines from 35 m to 1 km. The primary 
beam was ~ 2' at the line frequencies. The flux density scale 
was established by observing the standard primary calibrators 
3C286 and 3C48, and the uncertainty is expected to be less than 
~ 15%. Gain calibration was ensured by frequent observations 
of the compact quasar J0555+398. The quasar 3C84 was used 
for passband calibration. Pointing corrections were derived from 
X-band observations of nearby quasars and applied online. The 
data were edited and calibrated following standard tasks and pro- 
cedures of the Astronomical Image Processing System (AIPS). 
Imaging and deconvolution were performed using the IMAGR 
task, applying natural weighting to the visibilities. The resulting 
synthesised beam full widths at half-maximum are 2'.'33 x 2'.'27 
(position angle 92°) and 2'.'20 x 2'.'16 (position angle 90°) for the 
(1,1) and (2,2) lines, respectively. The noise level in each chan- 
nel is about 1.4 mJy beam"' and 1.5 mJy beam"' for the (1,1) 
and (2,2) data, respectively. 

Calibrated channel maps and spectra were analysed using the 
GILDAq3 software developed at IRAM and the Observatoire de 
Grenoble. 



3. Immediate results and derivation of physical 
parameters 

3.1. Integrated intensity maps, moments and averaged 
spectra ofNHj, (1,1) and (2,2) 

In the left panel of Fig. [1] we show the zero-th order moment 
map of NH3 (1,1) obtained towards 1, by integrating over the 
main group of hyperfine components. The shape is irregular with 
a main central filamentary structure extended > 80" in the NS 
direction, and a second clumpy filamentary-like structure to the 



' The VLA is operated by the National Radio Astronomy 
Observatory (NRAO), a facility of the National Science Foundation op- 
erated under cooperative agreement by Associated Universities, Inc. 

- The GILDAS software is available at |http://www.rram.fi7| 
IRAMFR/GILDAS 
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Fig. 1. Left panel: zero-th order moment map of the NH3 (1,1) line obtained with the VLA towards 1 (red contours). Contours start 
from the 6% of the maximum (0.039 Jy beam"' km s"'), and are in steps of 10% of the maximum (0.065 Jy beam"' km s"'). The 
integration velocity interval is -19.33 and -16.24 km s"', corresponding to the main group of hyperfine components. Contour levels 
start from the 3cr rms (~ 0.0012 Jy beam"' km s"'), and are in steps of 3 cr. The grey scale represents the 96 GHz continuum 
observed with the PdBI (first level is the 3cr rms = 0.004 Jy beam"' ; step = 3cr rms), and the compact 3 millimetre continuum cores 
detected in paper II are marked by crosses and labelled as Cl-a, Cl-b, C2 and C3. The two blue contours correspond to the 3cr level 
of the integrated emission of NaD^ (3-2) observed with the SMA (paper I). The arrows identify the emission peaks of the ammonia 
cores Al and A2 (this work). The star pinpoints the 3.6 cm continuum peak detected by Molinari et al. (2002). The dashed circle 
represents the VLA primary beam at ~ 22 GHz (~ 1 14"). The ellipse in the bottom left corner shows the synthesised beam of the 
NH3 (1,1) image. 

Right panel: Spectra of the NH3 (1,1) line towards the continuum cores Cl-a, Cl-b, C2 and C3, as well as the N2D^ condensations 
N and S, and the ammonia cores Al and A2 (see left panel) in synthesised beam temperature units {Tsb)- The spectra of Cl-a and 
Cl-b have been averaged over the 3cr rms contour level of the 284 GHz emission shown in Fig. 3 of paper II, because these are 
resolved in that image only. For C2 and C3, we used the 225 GHz 3cr rms contour level shown in Fig. 1 of paper I, because the 
225 GHz image is less affected by the flux filtering problem. 



NW having an arch-like shape (significance level from 6 to 9cr), 
oriented roughly in the NE-SW direction. The NH3 (1,1) in- 
tegrated emission traces an envelope that encompasses all the 
sources previously detected in the millimetre continuum and 
N2D"^ (Fig.[T] we adopt the same nomenclature as in paper II), 
in agreement with the relatively low critical density of the tran- 
sition (~ lO'* cm"-'). The largest integrated intensity is detected 
towards the millimetre continuum cores, i.e. in the direction of 
active star formation (see Sect. [U. 

Although it is difficult to identify separate cores given the 
complex emission shape, we point out two additional ammonia 
peaks labelled as Al and A2, seen also in N2H^ (see Sect. l3.2l i. 
Both peaks are undetected in the infrared, as we will show in 
Sect. 13.2.11 A2 is in between the continuum sources CI and C3 
in the central filament and located close to a marginal detection 
in the millimetre continuum (Fig. [TJ called C4 in paper II. C4 
is a marginal detection in the millimetre continuum, while A2 



is a well defined ammonia condensation, and it does not peak 
exactly at the position of C4. Therefore we cannot be sure that 
the two cores are associated. Note that we did not include C4 in 
any of the figures to make them clearer 

Fig.|2]shows the maps of the first- and second-order moments 
of the main group of hyperfine components of the (1,1) line. The 
first-order moment (panel a) does not show very steep velocity 
gradients. The second-order moment (panel b) clearly shows an 
enhancement of the velocity dispersion towards the continuum 
cores, i.e. close to active star formation. There is also evidence of 
line broadening towards the arch-like NW filament, around 30" 
north of the continuum source CI, with values up to 1 km s"' 
(corresponding to line widths of 2.35 km s"'). Both the velocity 
field and the velocity dispersion across 105345 will be examined 
more extensively in Sect. E] 

The NH3 (2,2) zero-th order moment map is shown in the 
left panel of Fig. [3] The emission is more compact than that of 
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Fig. 2. (a) first-order moment map (colour scale) of the NH3 (1,1) main group of hyperfine components. Same contour levels as in 
the left panel of Fig.[T]are shown. The synthesised beam of the NH3 (1,1) VLA observations is shown in the bottom-left corner. The 
symbols have the same meaning as in Fig.[Tl (b) same as panel (a) for the second-order moment. 
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Fig. 3. Same as Fig.[T]for the NH3 (2,2) transitions. In left panel, the NH3 (2,2) line emission has been integrated over the velocity 
interval -19.02 and -16.55 km s"' . Contour levels start from the 6% of the maximum (~ 0.015 Jy beam"' km s"'), and are in steps 
of 10%. 



the (1,1) line, and arises mostly from the mm-continuum cores. 
Significant integrated emission is detected towards the starless 
cores S, N, Al and A2. The arch-like filament to the NE is also 
clearly detected, but appears less extended than in the NH3 (1,1) 
map. For this line we do not show the first and second order 
moment maps which do not give any additional information with 
respect to that presented in Fig. |2] 



In the right panels of Figs. \T\ and [3] we show the spectra 
of NH3 (1,1) and (2,2), respectively, averaged over all identified 
cores. The contours used to extract the spectra are: the 3cr rms 
level of the 284 GHz continuum for Cl-a and Cl-b (paper II); 
the 3cr rms level of the 225 GHz continuum for C2 and C3 (pa- 
per I); the 3cr rms level of the NaD^ (3-2) integrated emission 
for N and S (paper I). We used the 3cr rms level of the NH3 (2,2) 
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emission to extract the spectra of Al and A2. Both Hnes are sys- 
tematically broader towards the continuum cores, showing once 
again that the other cores are more quiescent. 

The (1,1) lines were fitted by taking into account their hy 
perfine structure using METHOD NH3 of the CLASS packag^ 
For the (2,2) transition, the hyperfine fitting procedure did not 
give good results because the magnetic hyperfine components of 
the main central line were not resolved (expected maximum sep- 
aration of ~ 0.12 km s"'), and the inner satellites of the electric 
structure were not included in our observed bandwidth, so that 
we decided to fit the lines with single Gaussians. 

The results are listed in Table [1] For each core, small differ- 
ences are noticeable between Vpeak of the (1,1) and (2,2) lines 
beyond the indicated error margins (e.g., for CI, Cl-b, and N), 
probably because of the different fit method that does not take 
into account the hyperfine structure for the (2,2) line (see above). 
This could also explain the fact that Ad is found to be slightly 
but systematically smaller for the (1,1) than for the (2,2) line. 
However, in the protostellar cores this can have another origin: 
because the (2,2) transition needs higher temperatures to be ex- 
cited, it likely arises from gas closer to the forming protostars 
(as the integrated intensity maps in Figs. [T] and |3] suggest), and 
hence likely more turbulent. 

3.2. Comparison with otiier tracers 

3.2.1. Near- and mid-infrared continuum images 

The two ammonia filaments shown in Fig. [T] are located about 
~ 1' north-east of the near-infrared cluster corresponding to the 
IRAS source, as it can be seen from the 2MASS (/Tj-band) im- 
age (Fig.m. Two of the four IRAC images (band central wave- 
lengths 4.5 and 8 fim) of the region, taken with the Spitzer Space 
Telescope, are shown in Fig.|5] We do not show the 3.6 and 5.8 
fixn IRAC images because they are very similar to the 4.5 yum 
one. The star-forming cores Cl-a, Cl-b, C2 and C3 are all as- 
sociated with mid-infrared emission, while the starless cores N, 
S, Al and A2 lie in infrared-dark regions (see Figs. Eland |5)- 
In both the 2MASS and 4.5 fim image, we highlight the pres- 
ence of a point-like source located approximately 15" north of 
Cl-a (see yellow circle in Figs. |4]and|5]l, in between Cl-a and 
the north-western filament, which seems placed at the centre of 
a region almost devoided of ammonia emission. We will discuss 
better this source and its possible influence in the surrounding 
gas in Sect. 14.11 The nature of the (near- to far-)infrared emis- 
sion in 105345 will be discussed through deeper sensitivity and 
angular resolution images in a forthcoming paper (Fontani et al. 
in prep.). 

3.2.2. Higli density gas and outflow tracers 

The ammonia emission is in very good agreement with that of 
N2H"^ as can be seen in panel (a) of Fig. |6l where we superim- 
pose the integrated intensity of the main group of hyperfine com- 
ponents of N2H^ (1-0) (paper I) on the NH3 (1,1) emission map 
shown in Fig. [T] Because the primary beam of the N2H^ map 
is smaller that that of the NH3 maps, the N2H^ emission looks 



Table 1. Derived line parameters for the cores from the spectra 
integrated above the 3cr level. The uncertainties obtained from 
the fitting procedure are in parentheses. 



^ This method assumes that all the hyperfine components have 
the same excitation temperature and width, and that their separation 
is fixed to the laboratoiy value. The method also provides an es- 
timate of the optical depth of the line, based on the intensity ra- 
tio of the difi'erent hyperfine components. See the CLASS manual 
at http://i ramIr/IRAMFR/GILDAS/doc/html/class-html/class.htmr| for 
details. 
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3.3(0.6) 


-18.24(0.02) 


1.42(0.05) 


0.4(0.2) 


Cl-a 


8.5(1) 


-17.97(0.02) 


0.94(0.07) 
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Cl-b 


5.2(0.9) 
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1.3(0.1) 
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-17.78(0.02) 


0.87(0.06) 
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2.57(0.07) 
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3.11(0.03) 
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Notes. '"' Cols. 2-5 list the results of the CLASS hyperfine fitting pro- 
cedure: A X T„ = /[y,,(rex) - Jy(TgQ)], where / is the filling factor 
(assumed to be unity), JyiT^^) and yv(rBG) are the equivalent Rayleigh- 
Jeans excitation and background temperatures, respectively, t,„ is the 
opacity of the main group of hyperfine components; Vpeak = peak veloc- 
ity; Av = full width at half maximum corrected for hyperfine splitting; 
T,„ = opacity of the main group of hyperfine components; **' Cols. 2-5 
give the results of Gaussian fits: Area = total integrated intensity; Vpeak 



= peak velocity; Av = full width at half maximum; T^^^ 
sity in synthesised beam temperature units. 



peak 



: peak inten- 



less extended, but inside the PdBI primary beam (~ 48") the two 
tracers are almost coincident, although N2H^ (1-0) is expected 
to trace gas at slightly higher densities (wcrit ~ 10^ cm"-'). N2H^ 
(1-0) is clearly detected towards both Al and A2 indicating that 
these are real physical structures, even though in paper II we did 
not discuss them because they fall at the edge of the PdBI pri- 
mary beam. A more exhaustive comparison between NH3 and 
N2H+ will be presented in Sect. 15.11 

In panel (b) of Fig. |6] we plot the overlay between the NH3 
(1,1) emission and the widespread outflow detected in '^CO and 
reported in paper II. The internal part of the blue lobe of the out- 
flow coincides spatially with the elongated mid-infrared emis- 
sion detected west of Cl-b (compare panel (b) of Fig. |6]with 
Fig.|5]), suggesting that this emission could arise from the warm 
gas associated with the outflow. 

CARMA CS (2-1) observations of the region (Lee et 
al. 201 1 ) reveal seven emission peaks, which are plotted in panel 
(c) of Fig.|6l CS emission is detected towards CI, C3, Al, and 
close to the southern part of N. Cores C2, A2, S and N are not 
detected in CS. Because CS is well known to be highly depleted 
towards low-mass pre-stellar cores (e.g. Tafalla et al. 2002), the 
non detection in cores N and S could be due to molecular freeze- 
out, which would support the conclusions of paper I and II that N 
and S are likely in the pre-stellar phase. Core A2, not detected 
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Fig. 4. 2MASS /Tj-band image of the region associated with 
IRAS 05345+3157 (colour scale). The white contour represents 
the 6% level of the maximum (~ 0.039 Jy beam"' km s"') of the 
NH3 (1,1) zero-th order moment map. The yellow circle centred 
roughly at Dec. +32°00'20" pinpoints the location of the 4.5 yum 
source discussed in the text (Sect. 13.2.1b . The position error of 
the IRAS source is indicated by the grey ellipse. 



in CS, marginally detected in the milUmetre continuum (C4 in 
paper II), and clearly detected only in the non-depleted species 
NH3 and N2H^, might be another very young protocluster mem- 
ber where CS is depleted like in the pre-stellar core candidates. 
On the other hand, the non-detection of core C2 is puzzling since 
this condensation likely harbours the most evolved (and maybe 
most massive) object of the protocluster. 

3.3. Derivation of piiysicai parameters from ammonia 

3.3.1. Rotation temperature and NH3 column density 

From the averaged ammonia spectra shown in Figs.[r|and|3] fit- 
ted as explained in Sect. 13.1! we have derived several physical 
parameters. Rotation temperatures, Ti-oi, and column densities, 
A^(NH3), were calculated from the (2,2)-to-(l,l) line intensity 
ratio as described in Appendix A of Busquet et al. (2009). The 
method, based on Ho & Townes (1983), assumes that: (1) the 
inversion doublets constitute a two-level system; (2) the exci- 
tation temperature, Tex, and line widths at half maximum, AV, 
are the same for both transitions; (3) the physical conditions of 
the gas are homogeneous along the line-of-sight, so that if gradi- 
ents are present, the derived parameters are average values along 
the line-of-sight. In the derivation of the parameters, we also as- 
sumed a unity filling factor We stress that the T^oi derived using 
the results of Gaussian fits for the (2,2) transition (see Sect. 3.1) 
should differ only of 15% from the Ti-oi that we would derive by 
fitting the hyperfine structure. 

Trot and A^(NH3) are given in Table |2] As expected, the cold- 
est condensations are the pre-stellar core candidates N and S, 
with a Trot of 13 and ~ 16 K, respectively, while the continuum 
cores Cl-a, Cl-b, C2 and C3 have temperatures larger than 17 K. 
Because Cl-b harbours a hot-core (detected in paper II through 
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Fig. 5. Mid-infrared emission observed towards 105345 with 
the Spitzer Space Telescope. The two colour images show the 
Spitzer IRAC bands centred at 4.5 (top panel) and 8 jum (bottom 
panel). The white contour indicates the ammonia (1,1) emission 
and is the same as in Fig.|4] The yellow circle about 15" above 
Cl-a pinpoints the location of the 4.5 /^m source discussed in the 
text (as in Fig.|4|, and the other symbols indicate the location of 
the gaseous cores as in Fig. [1] The two grey contours represent 
the 3cr level of the N2D+ (3-2) integrated emission, and are the 
same as the blue contours in Fig. [1] In each panel. The VLA 
synthesised beam is shown in the bottom-left corner. 



CH3CN and OCS lines), the rotation temperature derived for this 
condensation is representative only of the lower-density enve- 
lope in which the hot-core is embedded. The source averaged 
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Fig. 6. Overlay between the NH3 (1,1) integrated emission (grey contour, same level as in Figs.|4]and|5]l, and the integrated emission 
of three previously observed tracers. They are, from left to right: (a) the main group of hyperfine components of N2H^ (1-0) observed 
with the PdBI (paper I, blue contours. First level and step = 0.02 Jy beam '); (b) the '^CO integrated emission over the red and blue 
non-Gaussian wings derived fro m SMA + NRAO merged data (same contours as in paper II); (c) the CS (2-1) line emission peaks 
detected with CARMA (Lee et al l201lL peak positions are taken from their Table 1). 

In each panel, black contours indicate the 3cr rms contour level of the N2D^ (3-2) emission, and crosses mark the peak position 
of the 3 mm continuum sources and of the ammonia cores as in Figs. [1] and [3] The ellipse in the bottom left corner and the dashed 
circle represent the synthesised beam and the primary beam, respectively, of the instrument used to observe the tracer superimposed 
on the ammonia emission: PdBI (left panel); SMA + NRAO (middle panel); CARMA (right panel). 



ammonia column densities are different among the cores, going 
from ~ 4 X lO''* cm^^ for Cl-b and C3 to ~ 1.3 x 10'^ cm^^ for 
C2. We will come back over this point in Sect. 15. II We stress that 
the uncertainties on the column densities given in Table|2] which 
are very small (less than 10%), do not take into account the error 
in the flux calibration, typically of about 10 - 15%. Therefore, 
an uncertainty of about 30% seems more realistic for A^(NH3). 

3.3.2. Core gas masses 

The gas mass can be estimated from the ammonia total column 
density from the relation: 



Table 2. Physical parameters from ammonia for the cores. 
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N(NH3) - 



41n2 
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A^(NH3) 
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where // - 2.33 is the mean molecular weight, ;«h is the hy- 
drogen mass, D is the linear diameter, A^(NH3) is the ammo- 
nia total column density and X is the [NH3/H2] abundance ra- 
tio. We adopt a uniform density as zero-th order assumption (i.e. 
A^(NH3) = n(NH3) x D) because we do not know the density 
profile of the condensations. 

To have an estimate of the NH3 abundance in our region, 
we can use the H2 column densities computed in paper II from 
the millimetre continuum emission (see Table 3 of paper II): we 
recomputed the masses from the optically thin millimeter con- 
tinuum emission using Eq. (1) of paper II, assuming the temper- 
atures derived from NH3 in this work (and a gas-to-dust ratio of 
100). For Cl-b, which harbours a hot-core, we have used both 
the ammonia temperature, representative of the lower-density 
envelope of the core, and that computed from the hot-core tracer 
CH3CN in paper II (i.e. ~ 200 K). This implicitely assumes 
coupling between gas and dust, which seems reasonable at the 
densities of our cores (> 10^ cm"-'). We derive X ^ \ x 10"^ 
and X ^ 0.85 x 10 '^ for cores CI and C2, respectively. In 
massive star-forming reg ions X is found to be of the order of 
10-^ to 10-^ (Harju et al. ll993L Wu et al. ,2006,) . We decided to 
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Notes. *"' Derived from interferometric observations of: millimetre con- 
tinuum (CI, Cl-a, Cl-b and C2, paper II), N2D+ (N and S, paper I), and 
ammonia (Al, A2 and C3 this work); **' from the 225 GHz contin- 
uum integrated emission (paper II), assuming that T^at equals the dust 
temperature and a gas-to-dust ratio of 100; ''' from the 284 GHz contin- 
uum integrated emission (paper II), assuming that Trot equals the dust 
temperature and a gas-to-dust ratio of 100; *'" from the 284 GHz con- 
tinuum integrated emission, assuming a hot core temperature of 200K 
(paper II) and a gas-to-dust ratio of 100. *''' core not well-defined in the 
ammonia maps. 



adopt X - 10"^. As in paper II, we do not consider the contin- 
uum sources close to or beyond the interferometer fields of view 
(namely C3 and C4). The various mass estimates are listed in 
Table|2l We point out that the dust masses derived from the mil- 
limetre continuum are affected by factors of uncertainty which 
are difficult to quantify, especially the dust mass opa city co effi- 
cient and the gas-to-dust ratio (see e.g. Beuther et al. l2002h . ex- 
pected to introduce factors of 2 or more in the uncertainty of the 
dust mass, which hence also affect both the H2 column density 
and the NH3 abundance. Therefore, both the dust mass estimates 
given in Table 2 and the derived ammonia abundances must be 
regarded with caution. 



Fontani et al.: The protocluster in IRAS 05345+3157 



Vlsr NH3(1.1) 



Vlsr NH3(2.2) 




20 



-20 



-40 




I -(-A2 



C3 If 



-17 



I 



-19 



20 10 -10 
Aa(arcsec) 



20 10 -10 
Aa(arcsec) 



Fig. 7. Map of the peak velocity of the NH3 (1,1) (left panel) 
and (2,2) (right panel) inversion transitions fitted as outlined in 
Sect. 13. II Crosses indicate the emission peaks of cores CI -a, CI - 
b, C2, C3, and A2 (see Fig.[T]). Solid contours mark the location 
of the N2D^ cores, (see Fig. [TJ- In both panels, the dotted con- 
tours range from -19.5 to -16.5 km s"', in steps of 0.3 km s"', 
and the ellipse in the bottom left comer shows the VLA synthe- 
sised beam. 



4. Maps of the physical parameters 

From the datacubes of both NH3 (1,1) and (2,2) we have ex- 
tracted spectra on a grid of spacing r.'6 x r.'6. These have been 
fitted as outlined in Sect. 13.11 From the fit results, we have de- 
rived maps of 105345 in several physical parameters, the discus- 
sion of which is the main subject of this section. 



4.1. Velocity field 

Maps of the peak velocity of NH3 (1,1) and (2,2) for the main 
filament are shown in Fig.|7] These maps are more accurate than 
that shown in panel a of Fig. [2] because they are based on the 
fitting method that takes into account the line hyperfine struc- 
ture. However, both maps agree quite well with the first-order 
moment map. The maps in Fig. |7]have been made interpolating 
the input data points on a regular grid. Because of this method, 
we recommend to regard with caution regions containing several 
points with faint emission, namely along the external edge of the 
emission and in the southern portion of the main filament. 

In Sect. l3.2.Tl we have remarked the presence of an infrared 
source detected in both the 2MASS and the Spitzer-IRAC im- 
ages (see Figs. |4] and |5]l, placed in between the main filament 
and the arch-like northern one at the centre of a region devoided 
of ammonia emission. Such a feature resembles the 'hole' pre- 
dicted by the model of Arce et al. (i2011i) for an expanding 
bubble (see their Fig. 5 and Sect. 3). In order to investigate 
this aspect in more detail, we computed position-velocity (p-v) 
plots along perpendicular cuts marked in the left panel of Fig. [8] 
centred on the infrared source. As seen in the figure, the two 
position-velocity diagrams show a U-like shape. In the Arce et 
al.'s model, an expanding shell would give a ring in the (p-v) 
plots (see their Fig. 5), wh ile we only see the redshifted half 
of it. However, Arce et al. (1201 ll) suggests that such a feature 
is consistent with a source driving the bubble slightly displaced 
towards the observer with respect to the surrounding molecular 
gas, so that we mainly see the gas which is moving away from 
us. This suggests that the infrared source is sweeping out the 
surrounding material through its winds/radiation, with an expan- 



sion velocity (difl^erence between the 'tip' of the U-like feature 
and the 'tails') of about 1 km s ' (see Fig. [8]l. In addition, the 
largest NH3 (1,1) line width in the region is found very close to 
this infrared source (see Figure |2] right panel), suggesting again 
the interaction of its winds/radiation with the dense surrounding 
gas. Thus, the more evolved members of the protocluster seem 
to be kinematically perturbing the starless cores, like condensa- 
tion N. This finding raises further questions about interactions in 
protoclusters, such as to what extent are the physical (tempera- 
ture and density), and chemical (C-bearing species, NH3 /N2H^ 
ratio, deuterium fractionation, etc.) properties of starless cores 
affected by more evolved cluster members. Thus, the 105345 re- 
gion harbours cores which are eminently suitable for assessing 
the influence of formed stars on starless condensations. 



4.2. Temperature and column density 



Following the procedure described in Sect. l3.3.Tl we have derived 
rotation temperature and total ammonia column density across 
105345. The results are shown in Fig.|9] We have excluded from 
the plots the arch-like filament, in which the (2,2) transition as 
well as the satellites of the (1,1) line are generally too faint to 
derive accurate fit results. In general we also excluded from the 
analysis the pixels for which the (2,2) line is not detected at a ~ 
4cr rms level, and we have used the same interpolation method as 
for Fig. |7] thus we recommend to regard with caution especially 
the values at the edge of the emission and towards N. In the left 
panel we see that Tiot is mostly in between 15 and 20 K along 
the main filament, especially close to CI (both CI -a and Cl-b), 
C2 and C3, while it is lower (less than 15 K) towards N and in 
the eastern part of S. 

Although the temperature map does not reveal big gradients, 
one can notice three temperature enhancements: (east of) Cl- 
b, north-east of C2, and east of A2. The temperature enhance- 
ment in Cl-b could be due to a combination of internal heating 
from the hot core (paper II) as well as heating by the passage 
of the outflow, which seems to be launched around this posi- 
tion. The enhancement north-east of C2 could be due to either 
external heating by C2 (although this is not too likely because 
there is no clear heating associated with C2) or to the red lobe of 
the widespread CO outflow (see panel (b) in Fig. |6]l. However, 
because this temperature enhancement is found at the edge of 
the emission, where the errors on Trot are of 4-5 K, we cannot 
be sure that this is real. The enhancement east of A2 could be 
possibly tracing an embedded source which would be related to 
the (faint) continuum emission detected at this position. Only 
high-sensitivity millimetre and infrared maps can help us to un- 
derstand the nature of this temperature enhancement. Note that 
the average temperature of the plateau that surrounds the tem- 
perature enhancements east of Cl-b and east of A2 is 15.8 K 
(with a Icr rms of 1 K), i.e. significantly lower than the maxima 
measured (20 and 22 K, respectively). Therefore, the tempera- 
ture enhancement proposed, except that north-east of C2, can be 
considered statistically significant. 

In general, the measured temperatures across 105345 are 
higher than those derived towards infrared-dark clouds (IRDCs), 
believed to be the earliest phases of massive star and cluster for- 
mation, in which T^ot is of the order of ~ 8 - 13 K (Ragan et 
al. 2011). They are also significantly higher than in clustered 
cores in low-mass star-forming regions such as Perseus (Schnee 
et al. 2009, Foster et al. 2009), where typical T^ot from ~ 9 to 
~ 15 K are measured. This indicates that the presence of already 
formed intermediate- and high-mass protostellar objects make 
the temperature of the dense gas systematically higher 
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Fig. 8. a) NH3 (1,1) integrated emission (white contours) overplotted on the 4.5 fim Spitzer IRAC emission. Contours start at 10% 
of the peak intensity, (39.03 Jy beam"' m/s), and increase in steps of 15%. b) Position-velocity (p-v) plot of the NH3 (1,1) emission 
along the cut at PA 70° and centred on (-5",13'.'4) with respect to the phase centre, with positive offsets increasing as shown by the 
arrow in panel a), c) idem as panel b) with PA 160°. Centre for both cuts is at position where the cuts intersect. In the p-v plots, 
contours start at -0.003 and increase in steps of 0.0015 Jy beam '. 
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Fig. 9. Map of rotation temperature (Trot, left panel) and NH3 to- 
tal column density (A^(NH3), right panel) derived from the NH3 
(2,2)-to-( 1,1) ratio as described in Sect. 13.3.11 In left panel, con- 
tour levels (dotted lines) range from 1 1 K to 20 K in steps of 3 K, 
while they range from 4xl0''*cm"^ to 3.2xl0'^cm"^ in steps of 
4xl0'^cm"^ in the right panel. Crosses and solid contours in- 
dicate the millimetre continuum cores and the N2D^ cores, re- 
spectively, as in Figs. [1] and |3] In both panels, the ellipse in the 
bottom left corner shows the VLA synthesised beam. 



The ammonia column density is on average of the order of 
~ lO'^cm"^ across the filament (right panel in Fig.|9]l, with the 
highest peak in between S and C2, and significant peaks towards 
C3, A2, south of N and east of C 1 -b. 



4.3. Non-thermal contribution to tine veiocity dispersion 

Assuming the gas to be in Local Thermodynamic Equilibrium 
(LTE) so that T^ot approximates the kinetic temperature, T^, we 
can compute the non-thermal contribution to the observed veloc- 
ity dispersion from the relation: cTnth = -J'^obs ~ "^tti where cTobs 

is the measured velocity dispersion of the core {- A(;/(81n2)'^^ 
for a Gaussian line profile with Ay the measured FWHM of the 
line) and cr^h is the thermal velocity disperion. For NH3, as- 
suming a Maxwellian velocity distribution, cr^h can be approx- 
imated by cr,h ^ 0.022 ^jT\(Tj km s"', with Tk =T,of In Fig.fTOl 
we plot both cTi^tit and the ratio cr„th/cr,j,, derived from the NH3 
(1,1) lines. As expected, the non-thermal contribution to the line 
widths is more pronounced close to the centre of the star forma- 
tion activity, especially towards C2 that likely harbours an Hii 
region, but also towards CI and C3. For all these cores 0-^^ ex- 
ceeds (Ttii of a factor larger than ~ 3.5. Towards N, S, A2 and 
the diffuse intracluster emission of the main filament, the ratio is 
mostly in between 2 and 3, suggesting that the gas is more qui- 
escent here but still with a dominant non-thermal contibution. In 
particular, we highlight the clear drop in the ratio o-nth/cth at the 
edge of core S, which resambles the sharp transition in the gas 
turbulence observed inside and out of the low-mass dense core 
Barnard 5 (Pineda et al. 2010). 

The plot is in agreement with the one made in paper II from 
the N2H^ (1-0) line widths (see left panel of their Fig. 11) which, 
however, was limited to the central portion of the main filament. 
The increase in the (Tnth close to YSOs is fou nd als o in low-mass 
star-forming regions (see e.g. Pineda et al. 1201 ll) . but the dif- 
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Fig. 10. Map of the non-thermal velocity dispersion (cr,ith, left 
panel) and ratio between non-thermal and thermal contribution 
to the observed velocity dispersion (cTnth/crth, right panel). In the 
left panel, the contour levels (dotted contours) range from 0.2 
to 0.6 km s"', in steps of 0.1 km s"', while in the right panel 
they range from 1 to 6, in steps of 1 . Crosses and solid contours 
indicate the millimetre continuum cores and the N2D^ cores, re- 
spectively, as in Figs.[3]and [T] 



feren ce here is more dramatic. Andre et al. (l2007b and Foster 
et al. (l2009l) found that in the low-mass clustered cores embed- 
ded in Ophiuchus and Perseus cr^tit is comparable to or smaller 
than the thermal broadening, respectively, regardless of whether 
the cores are isolated, clustered, starless or star-forming. On the 
other hand, in the few examples of protoclusters that contain 
more massive objects, line widths comparable to those measured 
in 105345 have been found (e.g. Palau et al. 12007 al Busquet et 
al.BoiOb, Filial etal. 2011). 

These findings suggest that on the spatial scales resolved by 
us, i.e. ~ 0.02 pc, the non-thermal motions are dominant in pro- 
toclusters that contain intermediate- and high-mass objects, and 
the dissipation of turbulence does not seem to occur at these sub- 
pc scales. High levels of turbulence on small spatial scales can 
be du e to th e overall high pressure of the environment (McKee 
& Tan 120021) and/or to large amounts of energy and momentum 
injected locally in the environment by the most massive proto- 
cluster members. The presence of the widespread '^CO outflow 
which covers a large portion of the dense gas and interacts with 
most of the condensations (see middle panel in Fig. [3]) suggests 
that the turbulence-injected scenario is very likely for 105345. 



5. Discussion 

5.1. NHi and /\/2H+ relative abundance: an evolutionary 
tracer? 

Unlike most of the C-bearing molecules, the nitrogen-bearing 
molecules NH3 and N2H^ do not suffer from depletion even at 
high-densities (i.e. in between ~ 10^ and ~ 10^ cm"^, e.g. Tafalla 
et al. l2002h and are therefore excellent tracers of dense gas . 
Studi es in b oth low- an d high-m ass star-forming regions (Hotzel 
et al. |2004 Palau et al. l2007ah suggest that the abundance ratio 
of the two species in dense cores is sensitive to evolution dur- 
ing the pre-stellar phase and the very beginning of the proto stel- 
lar phase, with the ratio NH3/N2H^ decreasing with time. Once 
CO is considerably desorbed from grain mantles by protostellar 



Table 3. NH3 -to-N2H^ column density ratio for the cores. In 
parentheses we give the error on the column densities. The un- 
certainty on the ratio A^(NH3)/A^(N2H+) is of the order of the 
60%. 



A'(N2H+) 



A'(NH3) 



(xlO'^cm-^) (xlO'^cm-^) 



A'(N2H+) 



CI protostar 1.3(0.4) 4(0.5) 30 

Cl-a protostar 0.8(0.2) 6.1(0.7) 76 

Cl-b protostar 1.2(0.3) 4.6(0.1) 38 

C2 protostar (Hii?) 2.8(0.4) 13(2) 47 

N starless (N2D+) 0.10(0.06) 6(1) 595 

S starless (N2D+) 0.23(0.05) 10(2) 427 

Al starless (no N2D+) 0.33(0.02) 2(0.3) 61 

A2 starless (no N2D+) 0.58(0.02) 6(1) 104 



feedback (either by outflow passage or by an increase in temper- 
ature), N2H^ is prog ressive ly destroyed and the ratio increases 
again (Busquet et al. l20IIh . However, the underlying physical 
and chemical reasons for the high values of NH3/N2H^ ratio 
during the pre-stellar phase and its subsequent decrease have not 
been investigated yet. 

To probe this intriguing chemical behaviour in the cores of 
105345, we have extracted spectra of N2H^ (1-Ojj and derived 
N2H^ column densities following the method described in pa- 
per I. We excluded C3 from the analysis because this is outside 
the primary beam of the N2H^ observations. In Table |3] we list 
the resulting column density ratio A^(NH3)/A^(N2H^). We find ra- 
tios from ~ 30 to ~ 80 in all the protostellar objects, ~ 60 and 
~ 100 for Al and A2, respectively, and ~ 400 and ~ 600 for 
the N2D^ condensations S and N, respectively. Because N and S 
are supposed to be pre-stellar, the relative abundance of NH3 to 
N2H^ seems greatly enhanced in pre-stellar cores (we point out 
that both Al and A2 were close to the edge of the primary beam 
of the N2D^, and were undetected, with a N2D^ column density 
upper limit of ~ 5 x 10' ' cm"^, corresponding to an upper limit 
on the N2D^-to-N2H^ column density ratio of ~ 0.01). 

That the NH3-to-N2H^ relative abundance in the gas phase 
is somehow enhanced in quiescent cores was firstly found 
in the intermediate-/high-mass protocluster IRAS 20293-1-3952 
by Palau et al. (12007a), in consistency with observations of 
lo w-mas s star-forming regions (Hotzel et al. 2004, Friesen et 
al. l2OI0h . In the low-mass regime, the NH3/N2H^ abundance ra- 
tio was also found to increase towards the centre of low-mass 
starless cores (e.g. Tafalla et al. 2002), which can be under- 
stood when freeze-out of species heavier than He becomes im- 
portant (e.g. Flower et al. 2006; Aikawa et al. 2005) . Depletion 
of heavy species like CO and other neutrals is also at the origin 
of the deuteration mechanism (see e.g. Millar et al.'l989') which 
could then explain the relation between high deuteration and 
high NH3/N2H^ abundance ratio. Given that no detailed expla- 
nation of the chemical processes driving the NH3/N2H^ increase 
can be found in the literature, we summarise the main passages 
in Appendix- A, where we also present the predictions of a sim- 
ple chemical model able to explain the observed NH3/N2H^ en- 
hancement in dense cores characterised by high values of molec- 
ular depletion. 



* from the PdBI observations published in paper I and using the same 
integration regions as those used to extract the ammonia spectra 
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5.2. Implications for the formation of (proto-)clusters 

An aspect with important implications for the formation of stel- 
lar clusters is how the observed separation between the members 
of a protocluster compares to typical fragmentation scale(s). A 
typical scale is the thermal Jeans length, which can be evalu- 
ated from t he temperature according to the relation (Stabler & 
Palla l2005h : 






-1/2 



(2) 



where c^ is the sound speed, andp and n are the mass density and 
volume number density, respectively. To compute a reasonable 
estimate of the Jeans length in i we consider an interval of 15- 
20 K for the temperature. For the density, the average density 
of the clump i n whic h 105345 is embedded is ~ 1 x 10* cm"^ 
(Fontani et al. l2006h . With these values, we obtain a thermal 
Jeans length Aj ~ 0.023 - 0.027 pc, corresponding to angular 
scales of ~ 2'.'5 . The spatial resolution of our observations is 
barely sufficient to resolve these scales. 

The average projected separation among the intensity peaks 
of the protocluster members within the central part of the main 
filament is 10" (see e.g. left panel in Fig. [TJ, larger than the ex- 
pected thermal Jeans length. The two sources separated more 
than the Jeans length from any other cluster member are A2 and 
C3. However, it is still unclear which is the nature of these con- 
densations and whether they are real protocluster members (es- 
pecially C3 which is quite isolated from the rest). If we consider 
only the cores that almost certainly belong to the protocluster, 
namely Cl-a, Cl-b, C2, N and S, their projected separations 
are roughly consistent with thermal fragmentation. On the other 
hand, the masses are of one order of magnitude higher than the 
expected Jeans mass (see Table |2]l, Ms - p/lj, estimated to be 
between ~ 0.1 and 0.18 Mq. 

Zhang et al. (2002) found a similar result in the IRDC 
G28. 34-1-0. 06, with core separations comparable to the Jeans 
length and core masses larger than the Jeans masses, and con- 
cluded that non-thermal motions such as turbulence and/or mag- 
netic fields can be dominant in the fragmentation of the IRDC. 
Similar findings have been found by Pillai et al. (,2011.) . who 
derived that the fragmentation in two IRDCs cannot be ex- 
plained by thermal fragmentation alone, but must include sig- 
nificant contributions from non-thermal motions. Because of the 
high non-thermal velocity dispersion measured (Sect. 4.3), a 
turbulent-driven fragmentation scenario or a significant contri- 
bution of magnetic fields can be applied to 105345 as well. 



6. Summary and conclusions 

We carried out VLA observations of the ammonia (1,1) and (2,2) 
inversion transitions towards the protocluster 105345. Ammonia 
is the best thermometer in dense gas and allows one to derive 
both temperature, column density and kinematics of the dense 
gas associated with the protocluster. The most direct results of 
this work are the following: 



- the ammonia emission is distributed in one main central fila- 
ment that encompasses all the previously detected high den- 
sity cores, and a secondary one to the NW of the protocluster 
centre; 

- the starless condensations have rotation temperatures of 
about 15 K, while the protostellar cores are hotter (Tiot in be- 
tween ~ 17 and 22 K). The temperature across 105345 is on 



average higher (roughly a factor ~ 1.3-1.5) than that mea- 
sured in both low-mass protoclusters and IRDCs using the 
same diagnostics. We suggest that this is due to the feedback 
of the already formed intermediate- to high-mass protostellar 
objects both inside and outside the protocluster; 

- the non-thermal contribution to the observed line widths is 
significantly higher than the thermal broadening both in the 
starless ((Tnth/cth ^ 2) and in the protostellar (cTnth/crth ^ 4) 
condensations. Specifically, the starless condensations are 
more turbulent than those observed in low-mass protoclus- 
ters (Ophiuchus and Perseus). Like the temperature enhance- 
ment, turbulence can be increased by the nearby massive pro- 
tostars and/or surrounding more evolved stellar clusters; 

- the NH3-to-N2H^ column density ratio seems greatly en- 
hanced in the pre-stellar core candidates N and S, which 
show NH3-to-N2H^ column density ratio a factor 10 larger 
than the others. Predictions of a chemical model show that 
this can be due to large freeze-out of molecular species heav- 
ier than He in dense cores. 

- By comparing the observed core masses and separations to 
the expected Jeans mass and Jeans length, we speculate that 
the turbulence played an important role in the initial frag- 
mentation of the parent clump, confirming similar evidence 
found in IRDCs. 
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Appendix A: NH3/N2H+ abundance ratio 

The NH3/N2H^ column density ratio has been found to signif- 
icantly increase in starless cores, in particular in those with the 
highest N2D^ to N2H+ abundance ratio (Sect. 15.1b . This trend 
is reminiscent of the NH3/N2H+ abundance increase towards the 
centre of low-mass starless cores (e.g. Tafalla et al. 2002), and 
can be understood in terms of freeze-out of species heavier than 
He (Flower et al. 2006; Aikawa et al. 2005). We present here the 
predictions of a chemical model, given that a detailed explana- 
tion of the chemical processes leading to the NH3/N2H^ increase 
is lacking in the literature. 

The freeze-out time scale (ffreeze) is inversely proportional 
to the density and becomes shorter than any dynamical time 
scale if the volume density (hh) is larger than ~ 10"^ cm""* 
(ffreeze - 10^/n/f yr; e.g. Bergin & Tafalla 2007). Therefore, in 
regions with «h larger than a few times 10"* cm""*, atomic and 
molecular species heavier than He becomes increasingly more 
depleted from the gas phase. In dark regions, this process is only 
counterbalanced by nonthermal desorption due to cosmi c-rays 
impulsive heating of dust grains (e.g. Hasegawa & Herbst [l993h 
and to photodesorption caused by UV photons produced by the 
interaction of cosmic rays with H2 molecules (Shen et al. l2004t 
Keto & Caselli 2010). Non-thermal photodesorption however 
does not depend on volume density, so freeze-out will eventu- 
ally win in the densest regions (riu > 10^ cm"^) of starless cores. 

Volatile elements such as H, D, and He are not affected 
by freeze-out. Actually, abundances of H^ (and its deuterated 
forms), He^ and H^ increase towards dense regions, where neu- 
tral species such as O and CO (their main destruction p artners ) 
are mostly frozen onto dust grains (e.g. Walmsley et al. l2004h . 



How does this affect the abundance of N2H^ and NH3? To an- 
swer this question, one needs to consider the main formation and 
destruction routes of these species. 

For this purpose, we used a reduced gas-phase chemical 
network extracted from KIDA0 and run models with different 
amounts of elemental freeze-out. The model includes H, He, O, 
C, N, S, Si, Fe and molecules up to 5 atoms in size. Initial con- 
ditions are based on standard low-metal elemental abundances 
(Table 1 of Caselli 2005), with all the elements in atomic form, 
except for hydrogen which is all in H2. The adopted H2 volume 
density, extinction and gas temperature are lO*" cm"-', 30 mag 
and 15 K, respectively. To study the effects of freeze-out, we 
varied the elemental abundance of all species heavier than He 
between a factor of 1 and 10"*. In Figure lA^ dop panel), the frac- 
tional abundace with respect to total H nuclei of key species and 
(bottom panel) the NH3/N2H^ abundance ratio (/?n) are plotted 
as a function of depletion factor (/d), together with the observed 
values from Table 3. Figure lA^ shows that 7?n first decreases 
with /d (when /d is between 1 and 7) and then constantly in- 
creases to values larger than 1000. The reason for this trend is 
due to the fact that N2H^ is directly linked to N2, given that its 
main formation route is 



N2 + Ht ^ N2H+ + H2 . 



(A-1) 



On the other hand, the formation of NH3 depends on N^, which 
can be produced by reactions of N-bearing species such as N2, 
CN, and NH2 with He+: 



N2 + He+ ^ N+ + N + He 



CN + He+ - 
NH2 + He+ 



C H- He H- N+ , 
■* He -H H2 + N+ 



(A-2) 
(A-3) 
(A-4) 



Reactions (IA-2b . ( IA-31 ) and (IA-41) are followed by successive hy- 
drogen abstraction reactions with H2 until all nitrogen valence 
electrons are used: 



N+ -> NH+ 



NH+ 



NHJ 



NHl 



(A-5) 



Once ammonium (NH^) is formed, it dissociatively recombines 
with electrons, to form ammonia. 

Reaction ( IA-2b dominates in non-depleted conditions, 
whereas reactions (IA-3b and (IA-4b dominate at large values of 
/d, when the formation rate of N2 drops. In fact, N2 forms via 
tw o (rela tively slow) neutral-neutral reactions (e.g. Hily-Blant et 
al. l2OT0l) : 



N + CH -> CN -I- H 

CN + N ^ N2 + C , 



(A-6) 

(A-7) 



so that an increase of /d affects the production rate of N2 more 
severely than the production rate of species formed via only one 
neutral-neutral reactions (e.g. CN, see reaction rA-6t . Moreover, 
as stated above, with the increase of /d, abundances of impor- 
tant destruction partners of N2 (H^ and He+) increases, further 
reducing the N2 abundance. This is evident in the upper panel 
of Figure lA-ll where N2 shows the steepest drop in abundance 
with /d. 

The N2H^ abundance increases between /d = 1 and 7, be- 
cause the drop in N2 is counterbalanced by a drop in CO, the 
main destruction partner at moderate values of /d (< 10): 



N2H+ + CO ^ HCO+ + N2 . 

^ Kinetic Database 

|http://kida.obs.u-bordeauxl.fr/| 



for 



(A-8) 
Astrochemistry; 
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At larger values of /d, dissociative recombination with electrons 
becomes the most important destruction route: 



NaH-" + e ^ N2 + H (90%) and N + NH (10%) 



(A-9) 



Reaction ( IA-81 ) reforms N2, which can produce again N2H+ via 
Reaction 



reaction (lA-lb . Reaction ( IA-9b forms N2 in 90% of the cases 
The other 10% of dissociative recombinations form NH, perma- 
nently removing N2 from the gas-phase. 

Destruction of NH3 happens via H^ (at large /d) and HaO^ 
(at low /d): 



NH3 
NH3 



m 



H30+ 



NH^ + H2 



nh; 



H2O 



(A-10) 
(A-11) 



However, unlike N2H'*', the destruction of NH3 at all values of 
/d reforms ammonium, the precursor of NH3. Therefore, the 
NH3 abundance drop follows the abundance of N^ precursors 
(see reactions IA-21 IA-31 and IA-4b 



From Figure lA-ll it is evi- 
dent that CN, NH2 and N (important for the formation of CN, 
see reaction lA-6b all maintain large abundances until /d ~ 1000, 
whereas N2H^ shows a constant drop for/o > 10. 
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Fig.A-1. Top panel: Fractional abundance of selected species 
as a function of depletion factor (see text for details). Note the 
sharp drop of N2 and the shallower drop of NH3 compared to 
N2H^. Bottom panel: NH3/N2H^ abundance ratio as a funciton 
of depletion factor. Horizontal lines are the observed values (see 
Table 3) for protostellar (continuous) and starless (dashed) cores. 
Note that although all protostellar cores are consistent with both 
low and large values of freeze-out, two of the starless cores can 
only be reproduce by values of depletion factor around 1000, 
similar to what measured in the core nucleus of LI 544 (Caselli 
et al. 1999). 
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